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ABSTRACT 

We have performed very large and high resolution cosmological hydrodynamic simula¬ 
tions in order to investigate detectability of nebular lines in the rest-frame ultraviolet 
(UV) to optical wavelength range from galaxies at z > 7. We find that the expected 
line fluxes are very well correlated with apparent UV magnitudes. The C iv 1549A 
and the C ill] 1909A lines of galaxies brighter than 26 AB are detectable with current 
facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C iv 
I549A, C III] I9O9A, [O ii] 3727A and [O ill] 4959/5007A can be good targets for the 
spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely 
Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space 
Telescope (JWST). Ha and H/3 lines are also expected to be detectable with these 
telescopes. We also predict detectability of the nebular lines for z > 10 galaxies that 
will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First 
Light And Reionization Explorer (FLARE) (11 ^ z ^ 15). We conclude that the G iv 
I549A, C III] I9O9A, [O III] 4959/5OO7A and lines even from z ^ 15 galaxies can 
be strong targets for TMT, ELT and JWST. We also find that the magnification by 
gravitational leasing is of great help to detect such high-z galaxies. According to our 
model, the C ill] 1909A line in z > 9 galaxy candidates is detectable even using the 
current facilities. 

Key words: cosmology - observations; galaxies - evolution; galaxies - formation; 
galaxies - high-redshift; galaxies - luminosity function, mass function; 


1 INTRODUCTION 

Understanding physical properties of the high-z galaxies be¬ 
yond z = 7 is very important, because such galaxies are 
not only building blocks of lower-z galaxies and local galax¬ 
ies but also the leading candidate of ionizing sources for 
the cosmic reionization. Such galaxies are much fainter ob¬ 
jects because they are very distant objects and less mas¬ 
sive systems. Thus, in order to observe these galaxies, ul¬ 
tra deep surveys are necessary. The Wide Field Camera 
(WFC3) installed on the Hubble Space Telescope (HST) 
in 2009 have detected many z > 7 galax y candidates us- 


ing the dropout 

Lyman break 

techniaue (lYan et al.l 

2011 

lEllis et al.ll2013l: 

McLure et al. 

I 2 OI 3 I: ISchenker et al.l 

2012 
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DunloD et al.ll2013l: iRobertson et al.ll2013l: lOno et al.ll2013l: 

Koekemoer et al. l2013l: 

Oesch et al. 20131: Bouwens et al.l 

2014al: Finkelstein et al. 

2014; Bouwens et al.l 2015ll. The 


combination of gravitational lensing by foreground galaxy 
clusters and the WFC3 camera have also dis covered even 
intrinsically fainter high- z galaxy candidates llOesch et al.l 
I2OI4I : llshigaki et akll^OlSl l . 


Many spectroscopic surveys for the Lya line of 
z > 7 galaxies have been attempted t o confirm spec¬ 


trosc o pic re d shift of these galaxie s (e.g., [PeAerkcTAY 


20121 

Shibuva et al.ll2012l: 

Finkelstein et al.ll2013l:lTreu et al. 

2013 

Schenker et al.ll2014 

:IVanzella et al.ll2014l;IOesch et al. 


201511 . Now, the highest confirmed redshift using Lya line 


is z = 8.68 JZitrin et al.ll2015bri which is the first and only 
example of a successful Lya line observation at z > 8. Other 
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observations so f ar have not succeeded in the Lya line de- 


tection at 2 > 9 (iBrammer et al. 

2 OI 3 I: iBunker et al.ll2013l: 

iGauak et al.ll2013l:lMatthee et al.l 

2 OI 4 J). Moreover, the spec- 


troscopic surveys for Lya line revealed that the fraction of 
high- 2 ; galaxi es with detectable Lya emission de c reases from 
z = 6 ('e.g.. IStark et al.ll201ll : lOno et al.ll2012l : iTiIvi et alJ 
I2OI4I '). These observational results imply that Lya attenu¬ 
ation by the neutral hydrogen in the intergalactic medium 
(IGM) is significant because the neutral fraction of the hy¬ 
drogen is high at such higher -2 before the completion of 
the cosmic reionization. The Lya line detection at 2 > 8 
may be very difficult even by future telescopes such as the 
James Webb Space Telescope (JWST), the European Ex¬ 
tremely Large Telescope (E-ELT), the Giant Magellan Tele¬ 
scope (GMT) and the Thirty Meter Telescope (TMT). Thus, 
it appears to be important to consider the detectability of 
alternative lines instead of the Lya line. Furthermore, the 
detection of high -2 metal lines helps to understand when and 
how the metal enrichment proceeds at very early Universe. 

The surveys for the far-infrared (FIR) lines using the 
Atacama Large Millimetre/submillimetre Array (ALMA) 
can be very powerful tool to confirm redshifts of very high- 
2 galaxies. The [C ll] 158/rm line is a well known bright¬ 
est line in the FIR region. The line has been detected in 


shifts (iMaiolino et al. 

200 a: 

lono et al.l I 2 OO 6 I 

l201ll:IVieira et al.ll2013 

:IWang et al.ll2013lL IG 


l2015lf . Interestingly, almost all survey have failed to de- 
tect [C ll] 158Mm line in high -2 Lyman a emitters 


jWyte^etalJj201^jKanekCTet^ 


Iota et al.ll2014l Watson et al.ll2015 L This fact disadvantages 


I 2 OI 3 I: lOuchi et all 12013: 


to explore high -2 [C ll] 158/im line because the fraction of 
young galaxies like LAEs increases with redshift. The [G 

II] 158/im is thought to be originated mainly from photo¬ 
dissociation regions (PDRs), the modelling of which is very 
complicated. This implies that the reliable theoretical pre¬ 
diction for t he [G ll] 158um line is very difficult. In our pre¬ 
vious work illnoue et al.l 120143 ). we concluded that the [O 

III] 88 /rm line from H ll regions is potentially the best line 
to measure redshifts of galaxies in the early Universe. The 
modelling of this line is simpler than that of the [G ll] 158/im 
line and the ALMA band 7 covers the line at 8.1 < 2 < 11.3. 

There are also various strong lines in the rest-frame UV 
to optical region. In the High-redshift Emission Line Survey 
(HiZELS) with narrow-band filters in the J, H and K bands, 
many Ha, [O ll] 3727A and [ 0 ill] 5007A emitters have been 
discovered at 0.4 < 2 < 2.2 (iGeach et al.ll2008l : ISobral et al.l 
l2009l . [ 2013 . 12015^ 1 . The spectroscopic surveys have also been 
performed and discovered many lines such as G ill] 1909A, 
[O ll] 3727 A, ]0 III] 5 OO 7 A and H a in the UV to the o ptical 
region at 2 < 3 (ISilyerman et al.||2014; iKashino et ahluOl^: 


Nakajimn^_&_OucM 2014 ■ Stark et al. l2014l : ISanders et al.l 
2 OI 5 I I. IStark et aiT ( 2015al ') detected G iii] 1909A emission 
in two high -2 galaxies at 2 = 6.029 and 7.213. The C 
III] I9O9A emission in unu sually luminous g alaxy at 2 = 
7.730 have been det ected bv IStark et al.l (l2016l L In addition, 
IStark et al.l (l2015bl ) discovered the G iv] 1549A and O iii] 
I665A emission lines in 2 = 7.045 galaxy confirmed by the 
Lya emission line. More rece ntly, a strong He ll 1640A was 
found in 2 = 6.6 LAEs ('CR7. ISobral et alJl2015 blL Int erest- 
ingly, the CR7 has no metal lines. Sobral et al.l ( ^15bll sug¬ 


gested the existence of metal-free stars so-called Pop HI stars 
and/or the formation of a direct collapse black hole. Now, 
the s ource of strong He ll 1640A in the CR7 is hotly debated. 

(e.g.j l^garwal et al.l2015l: Hartwig et al.ll2015l : IPacucci et al.l 


12015 !: Pallottini et al.r2(H5h. 


The UV to optical lines are very useful for studying 
the physical properties of galaxies such as the gas metal- 
licity, the radiation filed strength and the gas density (e.^ 


Nagao et al.ll2006l:lMannucci et ^l2010l : lNakaiima fc Oudiil 


20 14 Sanders et al. 2015l f . Nevertheless, there have been 


only a few line surveys at higher redshift (2 > 6 ) using the 
ground-based telescopes due to the observational difficulty. 
This is because the wavelength of these lines in the high -2 
Universe is redshifted into the IR region. The next genera¬ 
tion telescopes such as the JWST, the E-ELT, the GMT and 
the TMT can perform a deeper survey than the current ones 
in the IR regime. Thus, the metal lines can be good targets 
for these future telescopes and useful to diagnose the physi¬ 
cal properties of very high -2 galaxies (2 > 7). Unfortunately, 
only small number of bright galaxies at 2 > 7 appropriate for 
the follow-up survey can be discovered by these telescopes 
due to their narrow field-of-view (FOV) (~ 30 arcmin^). 
The Wide-Field Infrared Survey Telescope (WFIRST) and 
the First Light And Reionization Explorer (FLARE), which 
are also future telescopes, are the IR deep and wide-held 
imaging surveyors. The FOV of these telescopes are much 
larger than the TMT, the E-ELT, or the JWST. Many bright 
high -2 galaxies or very rare objects, which are good targets 
for spectroscopic follow-up with the TMT, the E-ELT, or 
the JWST, will be discovered. 

Many theoretic al studies have be e n made on very 
high -2 galaxies (e.g. , jFmlntor_etay| 2 p]T|jDa^aJ_etaJjl 20 lJ; 

I Shimizu et al. I [20141 : [ Glav et al. I I2OI5I: iFeng et al.ll2016h .1n 
these papers, they investigated the physical properties of 
star-forming galaxies at high -2 such as the halo mass, the 
stellar mass and the UV magnitude. They also revealed 
the evolution of the star-formation density, the stellar mass 
functions (SMFs), the UV luminosity function (LFs) and 
the metal enrichment in the IGM. There are some stud- 
ies for high -2 galaxie s with nebular emission consideration. 
ISchaerer et al.l (120091 ') studied nebular emission effects in the 
spectral energy distribution (SED). They showed that the 
estimated stellar mass, the age and the extinction strength 
chang e in th e cases of with or without the nebular emission. 
Ilnouel (1201 il l performed galaxy SED calculations in the UV- 
to-optical wavelength range with the nebular emission us - 
ing the photo-ionisation code cloudy dFerland et al.|[2013l l. 
In the paper, they discussed criteria of equivalent widths of 
some lines and broad-band co lours to select extrem ely metal- 
poor and metal-free galaxies. IWilkins et al.l (I2OI3II took into 
account the effect of the nebular emission on observational 
colour selections using a large cosmological hydrodynamical 
simulation. They claimed that the photometric redshift es¬ 
timation is strongly affected by the strength of the nebular 
emission. Although there are many studies for very high -2 
galaxies with the nebular emission, studies about the emis¬ 
sion line detectability in these high -2 galaxies are few. 

In this study, we examine whether some emission lines 
in the rest-frame UV to optical wavelength of galaxies at 2 = 
7-10 are detectable with the current and future telescopes. 
Next, we discuss the detectability of the lines at even 2 > 10 
in the future surveys. 


















































































































Detectability of nebular lines 3 


In section 2, we describe our numerical simulations and 
calibrate the parameters therein. We virtually observe the 
galaxies located in our simulation through a light-cone out¬ 
put and select simulated galaxies by applying a relevant 
magnitude limit. In section 3, we present some line fluxes in 
the rest-frame UV to optical wavelength of selected galaxies 
in our simulation. In section 4, we make predictions using 
the FLARE, the WFIRST, JWST and TMT telescopes. The 
final section is devoted to our conclusion. 

Throughout this paper, we adopt a ACDM cosmology 
with the matter density Hm = 0.3175, the cosmological con¬ 
stant Ha = 0.6825, the Hubble constant h = 0.6711 in the 
unit of Ho = 100 km s“^ Mpc“^ and the baryon density 
Hb = 0.04899. The matter d ensity fluctuations are nor - 


malised by setting as = 0.8344 (IPlanck CollaboratflonI 


All magnitudes are quoted in the AB system I Okd 


2014 1. 


199011 . 


The assumed initial mass function (IMF) in the observa¬ 
tional data and in our simulation is al ways the Chabr ier 
IMF with the mass range of 0.1-100 Mq llChabriejl2003ll . 



Figure 1. The stellar mass functions (SMF) at z = 7 and higher. 
The line styles and their corresponding redshifts are noted in 
the panel. For z = 7 and 8, we also plot observational dat a 
as the points with e rror bars taken f r om iGonzalez et all 1 I 2 OIIII . 
iGrazian et al.l ll2015^ and ISong et al.l 1 I 2 OI 5 I 1 . 


2 THEORETICAL MODEL 

In this section, we describe our cosmological hydrodynamic 
simulation, calculation of the spectral energy distribution 
(SED) of the simulated galaxies and the treatment of the 
dust attenuation for these galaxies. 


2.1 The Hydrodynamic Simulation 

We performed high-resolution numerical simulations with 
an updated version of the Tree-PM smoothed parti¬ 
cle hydrodynamics (SPH) code gadget -3 whi ch is a 
succe ssor of Tree-PM SPH code gadget -2 (ISpringell 
l2005h . We implemented relevant physical processes to 
galaxy formation such as the star formation, the 
supernova (SN) feed back and the chemical enrich- 


ment following pai 

Ders ( 

Okamoto. Nemmen & Bowerl l2008l: 

lOkamoto. & Frenk 

to 

0 

0 

CO 

: lOkamoto et al.ll2010lL In our ore- 


vious work I^Okamoto. Shimizu & Yoshidal 20141'). the radi- 


ation pressure and the AGN feedback were newly imple¬ 
mented. We consequently reproduced various observational 
quantities such as the stellar mass functions (SMFs), the 
cosmic star formation history, the galaxy downsizing, and 
the relation between stellar mass and metallicity from z = 4 
to z = 0. The details of these processes are found in the 
above references. 

Here, we briefly introduce settings of our simulations. 
We employ a total of 2 x 1280^ particles for dark matter and 
gas in a comoving volume of 50/i~^ Mpc cube. The mass of 
a dark matter particle is 4.44 x 10® ^“^Mq and that of a gas 
particle is initially 8.11 x 10® Mq, respectively. The soft¬ 
ening length for the gravitational force is set to be 2 h“^kpc 
in comoving unit. The SPH gas particles can spawn star par¬ 
ticles when they satisfies a set of standard criteria for star 
formation, and then, their mass is reduced. SPH particles 
around evolved star particles get the mass including metal 
elements due to SNe and the stellar mass loss in the AGB 
phase, resulting in the mass reduction of the star particles. 
In order to identify simulated galaxies, we run a friends-of- 
friends (FoF) group finder with a comoving linking length 
of 0.2 in units of the mean particle separation to identify 


groups of dark matter particles as haloes. Then, we find 
gravitationally bound groups of at least 32 total (dark mat¬ 
ter + SPH -I- star) particles as substructures (subhaloes) 
in each FoF group usi ng SUBFIND algorithm developed 
bv ISpringel et aH (1200iff . We regard substructures that con¬ 
tain at least 10 star particles as our simulated galaxies. We 
note that the lowest dark halo and stellar masses are around 
10® Mq and 10® Mq, respectively. 

In order to calibrate parameters in the code such as the 
SN feedback efhciency, we compare our model with the ob¬ 
served stellar mass function (SMF) of galaxies. Fig. [T] rep¬ 
resents our SMF at z = 7 to z = 15. We also plot the 
observational results (the points with error bars) at z = 7 
and z = 8 taken from Gonzalez et ^ 1 ([20i 3), IGrazian et al.l 
ll2015h and lSong et al.l |2015h . Our model reproduces the ob¬ 
served SMF within error-bars. At the massive end of SMFs 
(> 10*^® Mq), our model may underestimate the SMF be¬ 
cause our simulation box is not sufficiently large. On the 
other hand, our simulation tends to overproduce less mas¬ 
sive galaxies (< 10^ ® Mq) when comparing the observa¬ 
tional data at z ~ 7. This may indicate that the SFR of less 
massive galaxies in our simulati on is 2 times higher than 
real ones. In our previous work (IShimizu fc Inoual2013h . we 
argue that the shape of the SMF depends on the definition 
of the stellar mass, i.e., with or without the remnant mass as 
the stellar mass. The effect is important for low-z galaxies, 
but, it is negligible at the high-z discussed in this paper. 


2.2 SED Calculation of Simulated Galaxies 


After the identification of the simulated galaxies, we 
calculate their SEDs. The SED of each star particle, 
which has its own age, metallicity, and mass, is calcu- 
lated by using the popu l ation synthesis code pbgase 2 
jFioc fc Rocca-Volmerang3 IlflATl l. Then, we sum up the 
SEDs of the star particles composing a simulated galaxy to 
obtain the total intrinsic SED of the galaxy. We also adopt 
the nebular continuum and the Hydrogen recombination 
lines (Lya, Ho, H(l, and H 7 ) calculated by pbgase 2 with the 
Case B approximation. Moreover, the m etal lin e s emi tting 
from H II regions a re calculated based on llnou i (I 2 OIII') and 
llnoue et al.l ( 2014al i. in which expected line fluxes are cal- 
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culat ed by the photo-ionization code cloudy (iFerland et alJ 
l2013li . These line luminosity (Lune) are proportional to H/3 
luminosity {Luis), 

Tline = (1 — fesc)Cnne{Z)Lf!_p, (1) 


where /esc is the Lyman continuum escape fraction and 
Ciine(^) is the metallicity dep endent emiss ion efhciency for 
each lin e which is calculated in llnou idioni) and llnoue et alJ 
ll2014ah . respectively. Measurements of /esc of an order of 
0 . 01 - 0 .1 have been obt ained at 2 ~ 3 (e.g., Ilnoue et al.l 
I2OO5I : llwata et al.l l2009l ~l and /esc ~ 0 .2 at 2 > 5 has 
been favoured for c osmic reionization (Ilnoue et al. I I2OO6I: 
IShapley et al.| |2006l; ISiana et all |2015l) . However, we sim¬ 
ply assume /esc = 0 in this paper, yielding the maximum 
fluxes of nebular emission lines. It is worth describing ef¬ 
fect of a different choice of stellar population synthesis 
(SPS) model. One of important point for line emissivity cal¬ 
culation is the Lyman continuum photon production effi¬ 
ciency per nonionizing UV (~ 1500 A) because line emis¬ 
sivity strongly correlates wit h the efficiency. The difference 
between pegase 2 , starburst99 dLeithereret all 1999h. Maraston 
llMarastonll2005ll and bco3 (iBruznal &: Charlod 2003 ) is very 


small (e.g., Inoue et al.l l2014a : Wi]ldris_et_^ 


201(: ). How- 


ever, the efficiency of BPASS/binary (IStanwav et all 2015ll is 
abou t a factor of 2 higher than others fsee also lwilkins et al.l 
l2016h L Hence, uncertainty in the line flux estimation due to 
the choice of SPS models is also a factor of 2. 

After calculating the intrinsic SEDs of simulated galax¬ 
ies, then, we consider the dust attenuation for the continuum 
and the lines of them. The procedure of dust attenuation for 
the c ontinuum is same as our previ ous work dS himizu et al.l 
|2014 1. We adopt the Calzetti law (ICalzetti et al.l 2 OO 0 I I for 
the shape of the attenuation law. For the attenuation am¬ 
plitude, we calculate the escape probability of UV photons 
at 1500 A (//n?*). We apply the sandwich model for /{jy* 
dShimizu et a,l.l]20l3 : IXu fc Bu^Il995l l , 


/ uv ^ = ^(1 + + -(1 - e -"-^) 

Z Td 


( 2 ) 


where 5 is a parameter whose value is from 0 to 1 and Td 
is the UV (1500A) optical depth of the simulated galaxies. 
The parameter S is the fraction of the thickness of the central 
star+dust slab in the total thickness. For the case of <5 = 1, 
Eq.m corresponds to the well-mixed dusty slab geometry. 
On the other hand, it is the case with a central infinitely 
thin dusty sheet when <5 is zero. We calculate the optical 
depth Td using the following equation: 


3Sd 
4adS ’ 


(3) 


where ttd and s are a typical size and the material den¬ 
sity of dust grains, respectively. We set Od = 0.1 /rm and 
s = 2.5 g cm~^ motivated by SNe dust produ ction mod¬ 
els dTodini fc Ferarr"^ 1200 ll : iNozawa et al.ll2^003h . The dust 
surface density is dehned by the following equation: 


y Md 
= —2 > 

TTIT '^ 


(4) 


where Md and ra are the total dust mass and the effective ra¬ 
dius of the dust distribution in a galaxy, respectively. These 
values are directly obtained from our simulations as follows. 
The dust mass Md and are assumed to be proportional 
to the metal mass, Mmetai, and the half stellar mass radius. 


rhaif, of each simulated galaxy, respectively. Therefore, Md 
and rd are expressed as eMdMmetai and CrdPhaif, where eivid 
and Crd are the proportional constants of the dust mass and 
the effective radius, respectively. These two parameters can 
be reduced to one parameter. Thus, equation 0 is reduced 
to 


Ed 


6t 


Mmetai 


nr 


2 

half 


(5) 


where Ct = emd/e'^d is a global constant for all the galaxies 
in our simulation. We calibrate Ct and S so as to reproduce 
the observed UV luminosity function (LF) at 2 = 7, and we 
keep these values even at higher redshifts. In this study, we 
adopt Ct ~ 0.01 and S = 0.95, respectively. When we calcu¬ 
late band magnitudes, we also a pply the IGM absorption for 
the blue side of 1216 A following iMadaul lll995ll .li none et ^ 
ll2014bfl with an extrapolation for 2 > 7. The top panel of 
Fig.[2]represents the UV LFs from 2 = 7 to 2 = 10. Observa¬ 
tional results are also shown in the same panel. Our model 
reproduces the observations well up to 2 ~ 10. Moreover, we 
predict the UV LFs beyond 2 = 10 in the bottom panel of 

Fig.H 

Next, we consider the dust attenuation for various emis¬ 
sion lines. In this study, we apply the same attenuation 
law as the continuum to the emission lin es but with a 
different normalizati on (ICalzetti et al.ll 2000 l L According to 
ICalzetti et al.l ll 2000 h . the attenuation for emission lines is a 
factor of about 2 larger than that for the continuum at the 
same wavelength, altho ugh different facto rs were suggested 
by recent studies fe.g.. ICapak et al.ll2015ll . The line escape 
probability (/uv^) at 1500 A which is the normalization of 
the attenuation law for emission lines is 


rline 

/UV 


^Q-2.271og/5.-V"‘ 


(6) 


where /(jv* and 2.27 are the normalization of the attenu¬ 
ation law for the continuum (Eq. [5]) and the attenuation 
enhancement for emission lines, respectively. 

After calculation of the dust attenuation for the contin¬ 
uum and lines, we hnal ly calculate the dust thermal emission 
in t he same manne r as Shimizu. Yoshida fc Okamotd (l2012fl 
andlShimizu et al. l| 2014[l. One difference from our previous 
work ( Shimizu. Yoshida fc Okamotd l2012l : IShimizu et al.l 
l2014h is that not only the continuum photon energy but 
also the line photon energy absorbed by dust is converted 
into the dust luminosity. We assume that all the energy ab¬ 
sorbed by the dust is re-emitted in the IR emission. The 
monochromatic luminosity at frequency u is written 


= dTrMdK.HATd), 


(7) 


where Bi, and Td are the absorption coefficient, the 
Planck function and the dust grain temperature, respec¬ 
tively. We assume that 50 % of the total metal mass in a 
simulated galaxy is converted into the dust mass {Md) - The 
absorption coefficient in the FIR is well described by a 
power-law «:„ oc with /I = 1 ~ 2. In this s tudy, we ap¬ 
ply P = 1.7 based on obs ervational results dPunne et al.l 
120001 : IPunne fc EaleslI^Olfl . The typical dust temperature 
of our simulated galaxies is about 30 K. After these proce¬ 
dures, we obtain realistic SED of simulated galaxies from 
UV to IR/FIR wavelength. Fig. [3] represents an example 
SED from UV to IR of a simulated galaxy. As described our 
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Figure 2. The UV luminosity functions. The top panel represents 
the UV LFs from 2 = 7 to 2 : = 10. The bottom panel represents 
the UV LFs from 2 = 7to2 = 15. The line styles and their corre¬ 
sponding redshifts are noted in the panel. The points with error- 


llR.obertson et alJ 

2niCl: iBouwens et a 

.1 l2011allbl: iBradlev et alJ 

2 OI 2 I: Oesch et al, 

I 2 OI 2 I: iMcLure et al 

|201S 

: lOesch et alJ 20131: 

Schenker et alj l2013l: iBouwens et alJ 

2014a 

. I 2 OI 5 I: lAtek et alj 

2 OI 5 I: iMcLeod et alJ|2016|). 


previous work (jlnoue et al.ll201^ i. the [O ill] 88/rm line is 
the most prominent line. In the UV to optical range, many 
bright lines are also seen. We note that, in this study, the 
PDR/molecular cloud lines are not considered and the PAH 
features are not included in our SEDs, either. These are put 
into our future work. 


2.3 Basic Physical Properties of Simulated 
Galaxies 

Here, we explore whether our simulation is reasonable to 
predict line detectability. We check essential physical prop¬ 
erties of our simulated galaxies such as the stellar mass, the 
SFR, the dust attenuation and the metallicity which are 
necessary to calculate the line fluxes. In order to directly 
compare our model with the observations, we make galaxies 
distribution on a ligh t-cone from 2 = 6 to 2 = 17 the same as 
in our previous work llShimizu et al.l[2014l l. The field-of-view 
(FOV) in this study is (0.16 deg^) which is about 150 times 
wider than the HST FOV. We regard objects brighter than 
ffiso = 32 as our selected galaxies. Then, we divide them 
into four redshift subsets at 2 ~ 7 (6.5 <24 7.5), 8 (7.5 < 
2 < 8.5), 9 (8.5 < 2 < 9.5) and 10 (9.5 < 2 < 10.5). 


Figure 3. The spectral energy distributions of a simulated galaxy 
without the IGM attenuation. The solid, dashed, dot-dashed and 
dotted lines represent the SED of total components (stellar -|- neb¬ 
ular + dust), stellar component, nebular components (continuum 
+ line) and dust emission, respectively. 



Figure 4. Star formation rates as a function of the stellar mass. 
The large triangle, square, circle and star points are the median 
values of SFR in the stellar mass bins for 2 ~ 7, 8, 9 and 2 ~ 10, 
receptively. The boxes and error-bars ranges show 68 % percentile 
and 90 % percentile in the distributions. The observ ational data at 
2 ~ 6 also is shown as cross points with error-bars llSalmon et al.l 

I 2 OI 5 I) . 


2.3.1 Star Formation Rate 

The Star formation rate (SFR) of simulated galaxies is very 
important because the line luminosities of various elements 
are essentially proportional to the SFR. Fig. [4] shows the 
star formation rates as a function of the stellar mass. The 
large triangle, square, circle and star points are the median 
values of SFR in the stellar mass bins for 2 ~ 7, 8, 9 and 
2 10, receptively. The boxes and error-bars ranges show 

68 % percentile and 90 % percentile in the distributions. We 
also plot the observational data at 2 ^ 6 (cross points with 
error-bars) for comparison ll Salmon et al. 1 120151 ). The SFR 
strongly correlates with their stellar mass. Interestingly, the 
trend is similar to that of the star forming galaxies at lower- 
2 Universe known as the main seque nce distribution (e.g., 
iDaddi et al.ll2007l : iNoeske et al.l[2007l ) . We find that the dis¬ 
tribution of 2 ~ 7 galaxies is almost same as the observation 
at 2 ~ 6. Moreover, the evolution from 2 ~ 7 to 2 ~ 10 is 
very weak. This implies that line luminosities are also pro¬ 
portional to their stellar mass. 
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2.3.2 Dust Attenuation 

Fig. [5] represents the dust attenuation ^uv at 1500 A as a 
function of iifieo band magnitude. The large triangle, square, 
circle and star points are the median values of A\jy in i/ieo 
band magnitude bins for z ~ 7, 8, 9 and z ^ 10 , recep¬ 
tively. The boxes and error-bars ranges show 68 % percentile 
and 90 % percentile in the distributions. Clearly, there are a 
few galaxies strongly affected by the dust attenuation. The 
typical value of the dust attenuation correlates with the UV 
magnitude even though there is a large scatter. We note 
that the typic al value is smaller t han that we find in our 
previous work (IShimizu et al.l[2014l '). In the new code which 
we use in this study, the radiation pressure and the AGN 
feedback are newly included in addit ion to the SN feedback 
dOkamoto, Shimizu fc Yoshida|[2014l ~). Especially, the radia¬ 
tion pressure and the SN feedback effectively suppress the 
star formation in small mass galaxies than the case of the SN 
feedback only simulation. As a result, less dust attenuation 
is required to reproduce the UV LFs. 

Next, we calculate the UV slope P with the exactly same 
formulae as that in t he observations ijPunlop et al.l l2013l : 
iBouwens et al.ir2014bl : IWilkins et al.l[2016l l: 


P = 4.39 X (Ji25 — Hieo) — 2.0, ( 8 ) 

for z 7 galaxies, 

/3 = 8.98x (JHi4o-77i6o)-2.0, (9) 

for z ^ 8 galaxies, 

^ = 9.32x (JHi4o-77i6o)-2.0, (10) 

for z ^ 9 galaxies and 

P = l.lx {Hieo - [3.6]) - 2.0, (11) 


for 2 ~ 10 galaxies, respectively. J 125 and J77i40 are WFC3 
band filters equipped on the Hubble Space Telescope and 
[3.6] is the Spitzer IRAC bandl. Fig. [ 6 ] shows the slope 
P as a function of the absolute UV magnitude (without 
dust correction). The large triangle, square, circle and star 
points are the median values of P in the absolute UV mag¬ 
nitude bins for 2 ~ 7, 8 , 9 and 2 ~ 10, receptively. The 
boxes and error-bars ranges show 68 % percentile and 90 % 
percentile in the distributions. The la rge cross points are 
also shown as the observational data dDunlop et al.l l2013l : 
iBouwens et al.|[2014bl : I Wilkins et al.ll2016l l. As described our 
previous paper, the large scatter of the observational data 
may be due to the observational error and/or uncertainty. 
The UV brighter galaxies have larger P values. This is be¬ 
cause the UV brighter galaxies tend to be more evolved and 
have a smaller fraction of young stellar populations than 
the UV fainter ones. Moreover, such UV brighter galax¬ 
ies rapidly proceed metal enrichme nt. This trend can be 
seen in high -2 gala xy observations llBouwens et al.l l2014bl : 
IWilkins et al.l [2016l i. At 2 ~ 9, the P distribution of our 
model shows the large scatter and different trend from the 
other redshift ones. This is because Lyo break or line of 
some 2 ~ 9 (8.5 < 2 < 9.5) simulated galaxies enters in 
J77 i4o band. If we adopted the same P formula as 2 10 

for 2^9, we would get a similar P distribution for 2 ~ 9 to 
other redshifts. This suggests that some modifications are 
necessary for z ^ 9 P estimator. 



26 28 .30 32 

Apparent H^gQ magnitude [AB] 


Figure 5. The dust attenuation Ayv at 1500 A as a function of 
Hiqq band magnitude. The large triangle, square, circle and star 
points are the median values of Auv in Tlieo band magnitude 
bins for 2 ~ 7, 8, 9 and 2 ~ 10 , receptively. The boxes and 
error-bars ranges show 68 % percentile and 90 % percentile in the 
distributions. 

2.4 Metallicity 

In the fiducial line emission model of llnou iiEoH , the line 
emissivities are proportional to the metallicity in the case 
ol Z < 0.1 Zq where Z and Zq are the gas metallicity and 
the solar metallicity {Zq = 0.02). Then, when Z > 0.1 Zq, 
the dependence becomes opposite because of lower Lyman 
continuum (LyC) emissivity for higher metallicity. Thus, in 
order to detect the line emissions of very high -2 galaxies, 
the follow-up survey for ~ O.IZq systems is efficient. Fig. [7] 
represents the nebular metallicity of the simulated galaxies 
as a function of the stellar mass. The de finition of the metal - 
licity is the same as our previous work dShimizu et akllioillf 
in which the metallicity is dehned by the LyC luminosity 
weighted average metallicity called the nebular metallicity 
rather than the mass weighted average metallicity. The neb¬ 
ular metallicity of the simulated galaxies is proportional to 
their stellar mass. The metallicity of some galaxies reaches 
~ O.IZq even at 2 ~ 10. This means that the metal en¬ 
richment or the chemical evolution proceeds rapidly in the 
early Universe. Interestingly, the redshift dependence of the 
metallicity is very weak or almost no evolution. Thus, the 
each line emissivity at the redshift range varies little if their 
stellar mass is similar values. 


3 RESULTS AND DISCUSSION 

In this section, using our mock galaxies, we explore expected 
fluxes of some emission lines. We investigate the observed 
line fluxes of our simulated galaxies and discuss the de¬ 
tectability of these lines at 2 = 7-10. Then, we predict the 
detectability in 2 > 10 galaxies. 

3.1 Line Flux Prediction 

Here, we show some bright and important metal lines for the 
diagnostics of the gas metallicity. We also present hydrogen 
recombination lines such as Ha and Hp. The i723-index de¬ 
hned by ([O III] 4959A-B[0 hi] 5007A-b[0 ll] 3727A)/H^ is 
well known as a good indicator for gas metallicity. More¬ 
over, the N2 (log([N ii] 6584A/Ha)) and 03N2 (log(([0 iii] 
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Figure 6. The UV slope /3 as a function of the absolute UV magnitude. The large triangle, square, circle and star points are the 
median values of /3 in the absolute UV magnitude bins for 2 7, 8, 9 and 2 10, receptively. The boxes and err or-bars ranges show 
68 % percentile and 90 % percentile in the distributions. The large cross points are shown the observational results llDunlop et al.|[2013l : 
iBouwens et al.ir2014bl : Iwilkins et al.ll201^ 



* 2~/0 E 


10? ^qIO 

stellar mass [Mq] 

Figure 7. The nebular metallicity as a function of the stellar 
mass. The large triangle, square, circle and star points are the 
median values of the metallicity in the stellar mass bins for 2 ~ 
7, 8, 9 and 2 ~ 10, receptively. The boxes and error-bars ranges 
show 68 % percentile and 90 % percentile in the distributions. 
The dashed line is solar metallicity Zq = 0.02. 


5007A/H/3)/([N ii] 6584A/Hq) P indices are also used to es¬ 
timate the oxygen abundances. iNagao et all (l200(]l ~l suggests 
that the line ratio [Ne ill] 3869A/[0 ii] 3727A is also a use¬ 
ful metallicity indicator for high -2 galaxies, especially when 
the R23-index or other diagnostics involving [O ill] 5007A 
or [N ii] 6584A are not available. Fig. [ 8 ] represents the ob¬ 
served line fluxes of those important emission lines in the 
rest-frame UV to optical wavelength as a function of ob¬ 
served 77 i 60 band magnitude. We calculate all elements for 


the diagnostics (R23-index, N2 and 03N2) of the gas metal¬ 
licity. The small triangle, square, circle and star points are 
the expected line fluxes of the simulated galaxies at 2 7, 8 , 

9 and 10, respectively, while the large points with error-bars 
represent the average values and their standard deviations 
in 77 i6o band m agnitude bin s . We al so plot the o bserva ¬ 
tion al results by IStark et al.l (l2015al l. IStark et al.l (l2015bl l 
and lStark et al.l (I 2 OI 6 I I as the two cross and plus-sign points 
(the fainter point was demagnified by the quoted leasing 
magnification) in the Civ 1549A and the C ill] 1909A pan¬ 
els, which are very consistent with our prediction. 

We explore whether the existing facilities such as 
the VLT/X-Shooter and the Keck/MOSFIRE can detect 
the emission lines considered in this paper. The wave¬ 
length range of these instruments is about 1-2.5 /rm. 
Although, the detection limit is not deep enough (a 
fewxlO”^® [erg/s/cm^]) under a practical observational con¬ 
dition, the Civ 1549A and the C iii] 1909A lines of the 
brighter galaxies (< 26 mag) are detectable at 2 > 7 with¬ 
out a s ignificant magnificat ion by gravitational leasing. Re¬ 
cently, IZitrln et al.l ( 2015al l reported no C iii] 1909A line 
with ITieo 28 was detected under the condition of the de¬ 
tection limit 1.5 X 10“^® [erg/s/cm^j. This is also consistent 
with our result because typical line fluxes with Tfieo 28 
are 2 x 10 “^® [erg/s/cm^j which is below the observational 
detection limit. 

In order to detect these nebular lines, we need the future 
facilities with a wider wavelength coverage (the rest-frame 
UV to optical) and a stronger light-gathering power. The 
JWST and the TMT have potential to detect the nebular 











































8 Shimizu et al. 


line of very high -2 galaxies. We consider the detectability 
of the nebular lines using the JWST and the TMT. The 
JWST/NIRSPEC, the JWST/MIRI and the TMT/IRMS 
provide spectroscopy over the wavelength range of 0.6-5 /rm, 
5-28 fim and 0.8-2.5 fim, respectively. The point-source de¬ 
tection limits of the JWST/NIRSPEC, JWST/MIRI and 
TMT/IRMS, which depend on the observed wavelength are 
a fewxl 0 “^® [erg/s/cm^], a fewxl 0 “^® [erg/s/cm^] and a 
fewxlO”^'^ [erg/s/cm^], respectively, if the integration time 
is 1.0 X 10'* [sec] and the signal to noise ratio (S/N) is 5. 
We note that only C iv 1549A and the C ill] 1909A lines at 
2 = 7-10 can be available for the TMT/IRMS due to the 
wavelength coverage shorter than 2.5 ^m. Interestingly, the 
redshift evolution of the each line flux against the appar¬ 
ent Hieo band magnitude is not seen. This is because the 
SFR and metallicity evolutions of the simulated galaxies are 
weak as shown in Fig. [4] and Fig. [T] Moreover, the evolution 
of the dust attenuation is also weak as shown in Fig. [5l Al¬ 
most all these lines are detectable by these instruments if 
the band magnitude are brighter than about 28 mag. This 
means that we may detect some lines even at 2 > 10 galaxies 
if the galaxies are bright enough (< 28 mag) and we can con¬ 
firm the spectroscopic redshift using these lines. Moreover, 
we can study the gas metallicity of bright high -2 galaxies 
(< 28 mag) through the line ratio di agnostics. Finally, w e 
comment on Hell 1640A in the CR7 JSobral et al.l[2015bli . 
We also calculate the line and search the CR7 candidates 
around 2 = 6 . 6 . However, we can not find any CR7 like 
Hell 1640A bright galaxies in our simulation. This suggests 
that in order to reproduce such an object, we should con¬ 
sider the first sta r formation and/or intermediate black-hole 
formation (e.g., AEajwal_e^_alJ_ 20151: Hartwig et al.l l2015l : 


iPacucci et al.lT2015l : iPallottini et al.ll2015l i. This is beyond 

the scope of this paper. It will be our future work. 


3.2 Line Equivalent Width Prediction 


The equivalent width (EW) of lines is also very interesting to 
discuss detectability of lines in high -2 galaxies because the 
value strongly depends on the star formation history and the 
IMF. Fig. [9] represents the rest-frame EWs of various lines 
in the rest-frame UV to optical wavelength as a function of 
line fluxes. The large triangle, square, circle and star points 
are the median values of EWs in the observed line flux bin 
for 2 7, 8 , 9 and 2 10, receptively. The boxes and 

error-bars ranges show 68 % percentile and 90 % percentile 
in the distribu t ions. We also plot the obse rvati onal results 
Iw Stark et al.l (l2015af l. [Stark et ^ (l2015blf and lStark et al.l 
(l2016ll as the cross and plus-sign points in the Civ 1549A 
and the C ill] 1909A panels. The strong lines exceeding 
100 A can be seen in some lines such as ]Oiii] 4959/5007A, 
Hq and H/3. In the fainter line fluxes (less massive galaxies) 
regime (<'~ 10^° ]erg/s/cm^]), the EWs are proportional to 
their line fluxes. The EW represents the current SFR di¬ 
vided by an average of SFR over a certain time duration, 
and the line flux is almost proportional to the SFR. The 
metallicity effect on the line emissivity is secondary because 
the simulated galaxies have more or less a metallicity about 
0.1 solar (see Fig. [7|). If we choose galaxies emitting a 
weak line flux, they are in a low SFR phase, so that the EW 
is small. In the brighter line fluxes (massive galaxies) regime 
(>'~ 10^° ]erg/s/cm^]), EWs are almost constant (or very 


slowly decrease with increasing line fluxes). These galaxies 
are in an on-going star-formation episode (or starburst) and 
their EW (or SFR/ < SFR >) reaches an asymptotic value. 
Interestingly, only the Ha and H/3 lines show a decreasing 
EW as increasing the line flux rather than a constant. This 
is because stronger line flux galaxies tend to have a higher 
metallicity and the LyC photons production efficiency per 
nonionizing photon decreases with the metallicty increasing. 
On the other hand, for metal lines, the effect is cancelled 
out by an increase of the line emissivity due to increas¬ 
ing metallicity. Finally, we mention high er Ha EW galaxies 
tha n theoretica l expe ctation reported by I Shim et al.l (l201lf l 
and lSmit et al.l ll2016l l. In our simulation, there is no object 
which can reproduce such a very high EW. This implies that 
in order to reproduce such higher EW galaxies, we may need 
to include the metal-free star formation and/or a top-heavy 
IMF which we are not considered in this study. 


3.3 Detectability of the Lines beyond 2 10 

We have presented the detectability of multiple emission 
lines in the rest-frame UV to optical wavelength at 2 7- 

10. Here, we discuss the detectability of these lines at even 
higher- 2 , 2 = 11-15. The number of bright galaxies decreases 
with increasing redshift, thus, a very wide survey is neces¬ 
sary to discover sufficiently bright very high -2 galaxies. In 
this study, we feature high -2 galaxies which will be detected 
with the Wide-Field Infrared Survey Telescope (WFIRStJR 
and the First Light And Reionization Explorer (FLAREjj. 
The surveys with these facilities are potentially the most 
efficient one to detect very high -2 galaxies thanks to the 
wide survey area and the depth. We notice that the instru¬ 
ment design of these facilities has not been determined yet. 
In the survey strategy, they expect to detect many high -2 
galaxy candidates at 11 5/ 2 15. In order to compare our 

model with the future observations, we make a mock galaxy 
catalogue at 11 ^ 2 Si 15. We can identify a few galaxies 
at 2 15 even though the FOV of our simulation is very 

narrow (~ 0.15 deg^). 

Fig. [To] and [TT] represent the evolution of SFR and neb¬ 
ular metallicity, respectively. These results suggest that a 
weak redshift evolution of the line flux-magnitude relation 
even at 2 > 11 as lower -2 discussed in the previous section. 
Thus, some line emissions with K band magnitude < 26 mag 
can be detected even at 2 15 with the JWST/NIRSPEC 

and the JWST/MIRI. Next we estimate the line fluxes of 
these selected galaxies. 

Fig. [T5] represents the expected line fluxes in the rest- 
frame UV to optical wavelength as a function of K band 
magnitude. The point styles and their corresponding red- 
shifts are noted in the panel. The points with error-bars rep¬ 
resent the average values and the standard deviations in K 
band magnitude bins. The line flux-magnitude relations of 
different redshifts are almost the same as found in the cases 
of 2 ^ 10 galaxies as expected. We find that the C iv 1549A 
and the C ill] 1909A lines of simulated galaxies with K band 
magnitude < 28 mag are detectable even at 2 ~ 12 with the 
TMT/IRMS. On the other hand, in our simulation, there are 


^ http://wfirst.gsfc.nasa.gov/ 
^ http://mission.lam.fr/flare/ 
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Figure 8 . The observed line flux in the rest-frame UV to optical wavelength as a function of the apparent ffieo band magnitude. The 
large triangle, square, circle and star points are the median values of line fluxes in the i^ieo band magnitude bins for 2 : ~ 7, 8 , 9 and 
2 : 10, receptively. The boxes and error-b ars ranges show 68 % perce ntile and 90 % percentile in the distributions. The large cross and 

plus-sign points are observational results llStark et al.ll2015alla. [ 2 OI 6 I] . The solid, dot-dashed and dashed lines are the detection limits 
of the JWST/NIRSPEC (Aobs = 0.6-5 [p-m]), the JWST/MIRI (Aobs = 5-28 [/^m]) and the TMT/IRMS (Aobs = 0.8—2.5 [/^m]) at the 
wavelength of the line for a 2 : = 9 galaxy, respectively, if the integration time is 10“^ [sec] and the signal-to-noise ratio (S/N) is 5. 



Figure 10. Same as Fig. [T] but for 11 < 2 : < 15. 


no galaxies at 2 : ^ 15 that are bright enough. This implies 
that the survey area of our simulation (FOV^ 0.15 deg^) is 
not wide enough to have galaxies emitting line fluxes bright 
enough to be detectable with the JWST or the TMT. There¬ 
fore, a much wider survey area is imperative to discover line 
detectable galaxies beyond 2 11. We can detect some emis¬ 

sion lines using the JWST or the TMT if the bright galaxy 
candidates are discovered at any redshifts because there is 
no evolution in the relation between apparent UV magni¬ 
tude and line fluxes along the redshift. We also present the 
rest-frame EWs as a function of observed line fluxes in Fig. 



M Same trend as the 2 < 10 cases can be seen. There are 
very EW strong lines exceeding 100 A even 2 > 10. 

3.4 Nebular Emission Line Luminosity Functions 

We have not discussed yet how many line emitting galaxies 
can be detected by the future telescopes. This topic is very 
important for future plans. We show cumulative luminos¬ 
ity functions (LFs) of various emission lines at 7 ^ 2 ^ 15 
in Fig. 1141 Flere, we discuss only the C iv 1549A and the 
C III] lOOOA lines because the line have strong potential 
to be detected even at 2 > 10 with the TMT and the 
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Figure 9. The rest-frame EW of various lines in the rest-frame UV to optical wavelength as a function of observed line fluxes. The large 
triangle, square, circle and star points are the median values of EWs in the observed line flux bins for 2 ~ 7, 8, 9 and 2 : ~ 10, receptively. 
The boxes and error- bars ranges show 68 % perc entile and 90 % percentile in the distributions. The large cross and plus-sign points are 
observational results llStark et al.ll2015air3. [2016ll . 
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Figure 12. Same as Fig. [8] but for 11 < 2 : < 15. The solid, dot-dashed and dashed lines are the detection limit of the JWST/NIRSPEC, 
the JWST/MIRI and the TMT/IRMS at the wavelength of the line for a 2 = 12 galaxy, respectively, if the integration time is lO'^ [sec] 
and the signal-to-noise ratio (S/N) is 5. 
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Figure 13. Same as Fig. [9l but for 11 < ^ < 15. 
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JWST. Typical detection limits of TMT/IRMS (dashed 
line), JWST/NIRSPEC (solid line) and JWST/MIRI are 
also shown in Fig. [14] if the integration time is 10^ [sec] 
and the signal-to-noise ratio (S/N) is 5. The detection limit 
for TMT/IRMS is around 10“^® [erg/s/cm^]. Thus, one or 
more C iv 1549A and C ill] 1909A emitters at z < 12 can 
be detected by TMT/IRMS if the survey volume is about 
10® [Mpc®] which is equivalent to > 200 FOVs of the TMT 
and the JWST. On the other hand, the wavelength of the 
line at z > 12 is beyond the coverage of the instrument. 
The JWST/NIRSPEC can cover that wavelength range and 
the detection limit is around 10“^® [erg/s/cm^]. In this case, 
a larger survey volume is necessary to detect the emitters 
because there is no detectable galaxy in our small FOV 
(~ 0.15 deg^) simulation. We emphasise again that a wide 
imaging survey with the WFIRST and the FLARE is very 
useful to discover very higher-^ line emitters. 


3.5 Gravitationally-lensed z > 9 Galaxy 
Gandidates Found to Date 

Here, we mention detectability of galaxies magnified by 
gravitational leasing. Recently, many surveys with leasing 
magnification have been actively performed. Many faint 
galaxies have been discovered by the combination the 
HST/WFC3 and gravitational lensing. Thanks to the magni¬ 
fication by gravitational-lensing, some lensed galaxies which 
are candidates of 2 : > 9 objects are detectable by even 
the current facilities such as the VLT/X-Shooter and the 
Keck/MOSFIRE. 

Let us discuss the detectability of the C iii] 1909A 
emission line from 2 > 9 galaxy candidates found by sur¬ 


veys for galaxy cluster lensing fields like the Cluster Lensing 
And Supernova survey with Hubble (CLASH) and the Hub¬ 
ble Frontier Field (HFF) with instruments currently avail¬ 
able. The C III] I 9 O 9 A line is the second strongest emission 
line after Lya in the UV range and can be observed with 
ground-based 8-10 m class telescopes at 2 < 12. Table [T] is 
a summary of the candidates to be discussed here. First, 
we estimate the range of their real apparent Hieo magni¬ 
tudes by demagnifying with the quoted magnification factor 
jj, and its uncertainty. Next, we select the galaxies within 
the magnitude range and satisfying appropriate colour se¬ 
lection criteria from our simulation. Then, we obtain the 
magnification factor for each simulated galaxy so as to be 
observed as the observed Hieo magnitude. Finally, we mag¬ 
nify the C III] I 9 O 9 A emission line flux by this factor, and 
obtain the distribution of the expected C ill] 1909A line 
fluxes. Fig.[T5]shows the result for the five 2 > 9 candidates. 
The currently rea chable limit is 1 x 10“^® [erg/s/cm^] (e.g., 
IStark et al.ll20lJ l. The C iii] 1909A line from MACS1149JD 
is detectable with a high probability, 73 or 85% for the two 
quoted magnification factors. The detectability of the line 
from MACS0647JD1 is also relatively high (57%). The other 
three candidates have a detectability of < 5%. Therefore, 
MACS1149JD and MACS0647JD1 are the two best targets 
for the follow-up spectroscopy aiming to detect the C ill] 
I 9 O 9 A line even with the current facilities. In order to study 
higher -2 (2 > 7) line emitting galaxies, the combination of 
the GMT, the TMT or the JWST and gravitational lensing 
effect could be one of the most powerful ways. 
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Figure 14. The cumulative line luminosity functions at 7 < z < 15. The line styles and their corresponding redshifts are noted in 
the panel. Vertical lines represent typical detection limit of TMT/IRMS (dashed line), JWST/NIRSPEC (solid line) and JWST/MIRI 
(dot-dashed line) assuming the integration time is [sec] and the signal-to-noise ratio (S/N) is 5, respectively. 


Object 

■^ph 

Hieo [AB] 

h 

log 10(Fciii] 1909/[cgs]) 

P{< -18.0) [%] 

References 

MACS1149JD 

9.6 ± 0.2 

25.7 

± 

0.1 

14 c;+'^-2 / c. r-l-0.31 
-L"±- 0 _i 0 / ^-^-0.28 

-17.9 ± 0.1 

85/73 

ZI 2 , B14, H15 

MACS0647JD1 

10 . 20 - 11.47 

25.88 

± 

0.09 

8.4 

-18.0 ±0.1 

57 

C13, GL14, P15 

MACS1115JD 

Q 

^•^- 0.8 

26.2 

± 

0.2 

9 0+5.8 

-18.1 ± 0.1 

17 

B14 

MACS1720JD 

8.9t°j 

26.9 

± 

0.3 

= 0+4.7 

O.U_o 7 

— 18 4'*"^'^ 

0 

B14 

Abel2744JD 

9-8to.3 

27.37 

± 

0.16 

14 40'^^'^® 

Q6 

-18.6 ± 0.2 

0 

Z14, 115 


Tab le 1. A summary of gravit ationally-lensed 2 > 9 gal axy candidates. References: Z 12:[^eng et ^ l|201 jh B1 4: lBouwens et al.l' 2014a h 
H15 :lHu ang et al.| ||2015l)' Cl3: ICoe et al.l ll2013l) . GL14: [Gonzale2-Lopez et al.l ll20l4) . P15: |Pirzkal et al.l ||201^ , Z14: IZitrh^^iL ||2014 h 
115: llshigaki^^aLll20ld ^ 


4 CONCLUSION 

We have performed a large and high resolution cosmologi¬ 
cal hydrodynamic simulation to investigate the detectabil¬ 
ity of nebular lines in the rest-frame UV to optical wave- 
le ngth range at .g > 7. Our new simulation code based 


___ Qkamoto et al.l (|2010l) and lOkamoto, Shimizu fc Yoshid3 
(|20 i 4) can reproduce not only the stellar mass function, the 
downsizing, the mass-metallicity ratio but also the statisti¬ 
cal properties of star forming galaxies such as Lyman break 
galaxies (LBGs), Lyman a emitters and s u b-mm galaxies 
( Shimizu. Yoshida fc Okamotd I2OIII . l 2012 l : IShimi^ et alJ 
l2014h . First, we calibrate our model parameters so as to 
reproduce not only the stellar mass function at z 7 but 
also the UV luminosity function from z 7 to z 10. 
Then, we generate a light-cone output which extends from 
z = 6 to z = 17 using a number of simulation outputs. 
After making this light-cone output, we adopt the same 
colour selection criteria the same as in the observations with 


the HST/WFC3 camera dBouwens et al.ir2014al : [Oesch et al.l 
I2OI3II . Using the colour selected galaxies, we explore ex¬ 
pected emission line fluxes in the rest-frame UV to opti¬ 
cal wavelength such as C iv 1549A, C iii] 1909A, [O iii] 
5 OO 7 A and Ha, etc. These lines are very important not only 
to confirm spectroscopic redshift of the galaxies but also 
to perform the diagnostics of the gas metallicity. We find 
that the redshift evolution of the line flux as a function 
of the L7 i 60 magnitude is fairly weak because the metal¬ 
licity evolution is very weak as shown in Fig. [7] The C iv 
I 549 A and the C ill] 1909A lines from galaxies with the 
band magnitude < 26 AB at z > 7 are detectable even 
by the current telescopes suc h as the VLT/X-Sho oter and 
the Keck/MOSFIRE. Indeed, IStark et al.l l2015bli success¬ 
fully detected the lines from two galaxies at z ^ 7 and their 
detections are fully consistent with our prediction. Many 
important lines are detectable with the JWST/NIRSPEC, 
JWST/MIRI and the TMT/IRMS if the Hieo band magni¬ 
tude is brighter than 28 mag. Especially, the C iv 1549A, C 
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Figure 15. Probability distribution of the expected flux of the C 
III] I 9 O 9 A line for gravitationally-lensed z > 9 galaxy candidates 
found so far. For MACS1149JD, there are two estimates for the 
magnification factor and we show two cases (solid and dotted 
histograms). 


Ill] I 9 O 9 A, [O III] 4959 / 5 OO 7 A and H/3 lines are good targets 
for them. 

Finally, we predict the detectability of nebular lines for 
z ^ 11 galaxies satisfying K band magnitude < 30 mag. 
The redshift evolution of the line flux-K band magnitude 
relation is very weak even at 2 : ^ 11. We find that the C 
IV I 549 A and the C ill] 1909A lines are detectable by the 
JWST/NIRSPEC, and TMT/IRMS even at 2 - 12. How¬ 
ever, there are few galaxies at 2 ~ 15 to exceed the detection 
limit of the future telescopes in our small FOV simulation. 
This does not mean that we cannot detect emission lines 
from 2 ~ 15 galaxies. In fact, the JWST and the TMT can 
detect emission lines from galaxies with < 28 AB mag at 
any redshift which would be discovered, for example, by a 
> 100 deg^ survey planned with the future telescope such as 
the WFIRST and the FLARE. We can detect 2 ~ 15 galax¬ 
ies with a help of magnification by gravitational lensing even 
narrow FOVs. We predict that, according to our model, the 
C III] I 9 O 9 A line in 2 > 9 galaxy candidates (MACS1149JD 
and MACS0647JD1) is detectable using even the current fa¬ 
cilities such as the VLT/X-Shooter and the Keck/MOSFIRE 
with high probability. The survey with the combination of 
gravitational lensing and the GMT, the TMT, or the JWST 
is a good way to study higher -2 (2 > 7) line emitting galax¬ 
ies. 
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